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ABSTRACT 

We describe a numerical model constructed for the study of the emission of radiation from rela- 
tivistic plasma under conditions characteristic, e.g., to gamma-ray bursts (GRB's) and active galactic 
nuclei (AGN's). The model solves self consistently the kinetic equations for and photons, describ- 
ing cyclo-synchrotron emission, direct Compton and inverse Compton scattering, pair production and 
annihilation, including the evolution of high energy electromagnetic cascades. The code allows calcu- 
lations over a wide range of particle energies, spanning more than 15 orders of magnitude in energy 
and time scales. Our unique algorithm, which enables to follow the particle distributions over a wide 
energy range, allows to accurately derive spectra at high energies, > 100 TeV . We present the kinetic 
equations that are being solved, detailed description of the equations describing the various physical 
processes, the solution method, and several examples of numerical results. Excellent agreement with 
analytical results of the synchrotron-SSC model is found for parameter space regions in which this 
approximation is valid, and several examples are presented of calculations for parameter space regions 
where analytic results are not available. 

Subject headings: galaxies: active — gamma rays: bursts — gamma rays:theory — methods: numerical 
— plasmas — radiation mechanism: Non thermal 



1. INTRODUCTION 

In the standard fireball scenario of gamma-ray bursts 
(GRB's), the observable effects are due to the dissi- 
pation of kinetic e nergy i n a highly r elativ is tic fireball 
(see, e.g., iPiran I i200Cjk iMeszaros I l|2002fl : iWaxman I 
(2003) for reviews). Synchrotron emission and inverse- 
Compton emission by accelerated electrons are the 
main radiative processes. Electrons accelerated in 
the internal shock waves within the expanding fire- 
ball produce the prompt 7-ray emission, while elec- 
trons accelerated in the external shock wave driven 
by the fireball into the surrounding medium pro- 
duce the afterglow emission, from the X to the radio 
band s. " •' ' : ' r>1 ' x ' ' ' ' " 
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Paczvhski & Rhoads 1993; Meszaros & Rccs 
etri I Il997at iSari. Piran fc Naravan I I199F 



Gruzinov fe Waxman 1 11999') . 



n g in g eneral agreem ent w i th obse rvations 
ifiMlPreece et al. llT998l iFrontera et al. I 
2000; Meszaros 2002) both theoretical arguments and 



Whi le being in general agreement w i th observations 
Bande^^ 

TUT 

observational evidence suggest that the optically thin 
synchrotron - synchrotron self Compton (SSC) emission 
model is not complete in explaining neither the prompt 
nor the afterglow emission. Additional physical pro- 
cesses can significantly modify the SSC spectrum. First, 
over a wide range of model parameters, a large number 
of e ± pairs are produced in internal collisions, due to 
annihilation of high energy photons. Second, relativis- 
tic pairs cool rapidly to mildly-relativistic energy, where 
their energy distribution is determined by a balance be- 
tween emission and absorption of radiation. The emer- 
gent spectrum, which is affected by scattering off the pair 
population, depends strongly on the pair energy distri- 
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bution, and in particular on the "effective temperature" 
which characterizes the low-end of the energy distribu- 
tion. Third, proton and electron acceleration to high 
energies initiates rapid electro-magnetic cascades. It is 
necessary to follow the evolution of a high energy, non- 
linear cascade in order to accurately derive the spectrum. 
And last, the plasma is not in steady state, and the par- 
ticle distributions are continuously evolving. 

On the observational side, we note that hard spec- 
tra, vF v cx v a with a > 4/3 at low ener- 
gies, < 300 keV , were observ e d at early times in 
some GRB's llPreece et al. Ml998HFrontera et al. 11200(1 
iGhirlanda. Celotti fc Ghisellini 1 2003). These hard spec- 
tra is inconsistent with the optically thin synchrotron- 
SSC model predictions. A similar conclusion was ob- 
tained by a comparison of the high energy and low en- 
ergy spectral in dices during the prom pt emission phase 
of 150 GRB's l)Preece et al. I 120021) . An additional 
high energy (> 1 MeV ) spectral component inconsis- 
te nt with the synchrotron model prediction was reported 
bvlGonzalez et al 1 120031) . Finally, a recent analysis by 
iBaring fc Brabvl |2004) showed difficulty in explaining 
the high energy component of GRB's early emission by 
the SSC model. 

The above mentioned difficulties, raise the need for a 
model that can better describe emission under conditions 
characterizing GRB's. However, a numeric calculation of 
GRB spectra that takes into consideration creation and 
annihilation of pairs is complicated. The evolution of 
electromagnetic cascades initiated by the annihilation of 
high energy photons occurs on a very short time scale. 
On the other extreme, evolution of the low-energy, mildly 
relativistic pairs, which is governed by synchrotron self- 
absorption, direct and inverse Compton emission takes 
much longer time. The large difference in characteris- 
tic time scales poses a challenge to numeric calculations. 
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Another challenge to numerical modeling is due to the 
fact that at mildly relativistic energies the usual syn- 
chrotron emission and IC scattering approximations are 
not valid, and a precise cyclo-synchrotron emission, di- 
rect Compton and inverse Compton scattering calcula- 
tions are required. 

Two approaches have been employed so far in treat- 
ing this problem. The first is the Monte-Carlo ap- 
proach, where individual particles are followed as they 
undergo interactions inside the plasma. This scheme 
typically suffers from relatively poor photon statistics at 
high energies, and does not lend itself to time-dependent 
calculations. Work done so far using this approach 
ijPilla fc Loeb lll99S|) was limited to parameter space re- 
gion where the creation of pairs has only a minor effect 
on the resulting spectrum. The second approach involves 
solving the relevant kinetic equations. Following the time 
evolution of the system using this method is straightfor- 
ward, and photon statistics at high energies is not an 
issue. However, the above mentioned complications lim- 
ited ti£jic£ura£v_j3ft^^ constructed so 
far ijPanaitescu fe Meszaros lfl098f) above ~ 1 GeV . 

Note that this method was extensively used in the 
past in the study of active galactic nuclei (AGN) 
plasma llZdzia rski fc Lightma nT 119851 [Fabian et al. I 
Il986t iLiehtman fc Zdziarski lll987tlCoppi Ill992|h Using 
the numerical models new results were obtained, such 
as the effective pair temperature and the complex pat- 
tern o f the spectral indices in the X -ray (2 — 10 keV ) 
range ijLightman fc Zdziarski lll987lK that were not ob- 
tained by previous analytic calculations. Non of these 
models, however, considered the evolution of high en- 
ergy electro-magnetic cascades, expected to be relevant 
for both GRB's and AGN's. In addition, the treatment 
of photon emission in the presence of magnetic field was 
not complete, since particles are expected to accumulate 
at low energies (7 ~ 1), where the synchrotron emission 
approximation used does not hold, and exact treatment 
of cyclo-synchrotron emission is required. 

Pair cascade e v olution was first studied by 
IBonometto fc Rees I (j 1 9 7 1^1 . Small angle cascade 
showers in anisotr o pic ra diat ion field were treated by 
IBurns fc Lovelace I lIlQSl. ICuilhert. Fabian fc Rees I 
(1983) and iSvenssoiil l)1987t) have generalized the 
treatment of cascade evolution, showing that it may 
have a significant effect on the high energy spectrum. 
It is therefore necessary to incorporate the cascade 
calculation in order to accurately derive the high energy 
emission spectrum. 

We have constructed a numerical model that over- 
comes the numerical challenges. Applying this model to 
GRB plasmas, we have o btained several new results . For 
example, we have shown llPe'er fc Waxman"ll2004b|) that 
emission peaks at ~ 1 MeV for t± ~ 10 — 10 2 , where 
t± is the optical depth to scattering by pairs, and that 
peak energy at 3> MeV cannot be obtained for GRB 
luminosity L ~ 10 52 ergss _1 . We showed that for large 
compactness, I' > 100, the spectral slope below 1 MeV is 
steep, e 2 n p h{e) oc e Q with 0.5 < a < 1 and shows a sharp 
cutoff at 10 MeV . We also showed ijPe'er fc Waxmanl 
I2004cj) that observations of the early afterglow emission 
at 1 GeV — 1 TeV is informative about two of the most 
poorly determined parameters of the fireball model: the 
ambient matter density, and the fraction of thermal en- 



ergy carried by the magnetic field, es. 

We present in this paper our numerical model. In SjSJwc 
describe the basic model assumptions. We then present 
the kinetic equations that are being solved, and detailed 
description of the numerical treatment of various physi- 
cal processes. Our numerical integration approach is de- 
scribed in 33 We present the general approach of treat- 
ing this complicated problem, and the various integration 
techniques used. In <2]we give examples of numerical re- 
sults, relevant to the prompt emission phase of GRB's, 
and compare them to approximate analytic results. We 
summarize in §S| the main features of our method, and 
discuss its usefulness for the ongoing research of GRB's 
and AGN's. 

2. MODEL ASSUMPTIONS AND PHYSICAL PROCESSES 

We consider a uniform plasma, composed of pro- 
tons, electrons, positrons and photons, and permeated 
by a time independent magnetic field. The particle 
and photon distributions are assumed homogeneous and 
isotropic. Considering the physical phenomenon of, e.g., 
GRB as a motivation, these assumptions are equivalent 
to the assumption that the calculations are carried out in 
the comoving frame (see jJH below). We assume the exis- 
tence of a dissipation process, (e.g., collisionless shock 
waves) which produces energetic particles at constant 
rates, Q(^) and 5(7) electrons and protons respectively 
per unit time per unit volume per unit Lorentz factor, 
7. Since the details of the acceleration process are not 
yet known, we do not specify here the functions Q(~f) 
and S("j). These functions are specified when treating a 
particular problem (see . Motivated by the GRB fire- 
ball model scenario, in which the internal shock waves 
cross the colliding shells at relativistic speeds, we as- 
sume that the dissipation process occurs on a character- 
istic time scale which is equal to the light crossing time, 
tdyn ~ R/c, where R is a characteristic length scale of 
the plasma. 

The population of electrons, positrons and photons is 
affected by synchrotron emission, synchrotron self ab- 
sorption, Compton scattering, pair production and pair 
annihilation, that occur simultaneously during the dy- 
namical time, tdyn- In this version of the code, protons 
are assumed to interact via photo-meson interactions 
only, producing pions that decay into energetic photons 
and positrons. Coulomb scattering is not considered, be- 
cause, as we show in fflXfr. it is insignificant in calcu- 
lating the spectra u nder conditions that a re of interest 
to us. As noted bv iCoppi fc Blandford I l)1990|) . e — e 
bremsstrahlung is also insignificant under the same con- 
ditions, and is therefore not included in the calculations. 

We assume no photon escape during the dynamical 
time, tdym and instantaneous photon release at the end 
of the dynamical time. This approximation is justified 
since the dynamical time is equal to the light crossing 
time. If tt, the optical depth to Thompson scatter- 
ing by electrons or by the created pairs at the end of 
the dynamical time is larger than 1, the 'instantaneous 
release' approximation is not valid. Since the plasma 
is assumed to be heated to relativistic energy density, 
we assume in this case that the dissipation phase is fol- 
lowed by a relativistic expansion phase, during which 
the optical depth decreases. The evolution of particle 
and photon distributions is followed during the expan- 
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sion phase until the optical depth for Thompson scatter- 
ing, tt < 1. A detailed description of this calculation is 
given in lPe'er fc Waxmanl l)2004rj|) . 

Let n e - (7, i), n e + (7, t) and n p (7, t) be the number den- 
sity per unit Lorentz factor, 7, per unit volume of elec- 
trons, positrons and protons, and n p h(s,t) be the num- 
ber density per unit energy per unit volume of photons. 
The time derivatives of the electron, positron, proton and 
photon number densities are given by 



dn e -(n,t) d 



n B -(7»*) (Psh) + Pc{-r,t)) 



+H{n,t)M 



2 d ( n e -(7,t) 



+ 



dn e p(j,t) dn eA (-f,t) 



at 



Of 



(i) 



dn e + (7, t) 

at 



d_ 



n e+ ( 7 ,t) {Ps(l) + Pc(rt,t)) 



d ( n e +(7,t) 



(?7 V (3j 2 



dn eP (j,t) dn eA (j,t) 



dn p (j,t) 



dt 



<:)t 



at 



97 



[n p (7,t)P 7r ( 7 ,t)] 



(3) 



dn ph (e,t) 



— ^ — - R s (s,t) + R c (e,t) - R P (s,t) + R A (e,t) 
+i?7r(c,t) - cn ph (e,t)a(e,t). 

(4) 

Here, the terms in the parenthesis on the right hand 
side of equations 1 1121 give the change of population due 
to synchrotron emission and Compton scattering, Ps(j) 
and Pc{"j,t) are the synchrotron and Compton emitted 
power. The third term in the parenthesis represents en- 
ergy gain by synchrotron self absorption, with H{^,t) 
defined below, see equation 03 The last two terms in 
equation^ dn e p(j, t)/dt and dn e A{j, t)/dt are the rates 
of pair creation and pair annihilation per unit volume. 
The term (3^(7, t) in equation [21 represents positron cre- 
ation by the decay of tt + . The pions are produced by 
photo-meson interactions of low energy photons with en- 
ergetic protons. In the proton equation, Pn(j,t) is the 
rate of proton energy transfer to pions. In the photon 
equation, Rs(e,t), Rc(e,t), Rp(e,t) and RA{£,t) are the 
rate of production and annihilation of photons due to 
synchrotron emission, Compton scattering, pair produc- 
tion and pair annihilation, and i? T (e, t) is the production 
rate of photons due to the decay of energetic pions. The 
last term represents photons reabsorption, where a(e, t) 
is the self absorption coefficient. 

In the rest of this section, we present detailed descrip- 
tion of the terms in equations 

2.1. Synchrotron and synchrotron self absorption 
emission terms 

The term H (7, t) in equation ^ describes the heating 
of the electrons and their diffusion in energy due to syn- 



chrotron self absorption. It is given by 



duj 



(5) 



Csee IMcCravl 11963 iGinzburg fc Svrovatskii I IT963 
IGhisellini. Guilbert fc Svensson II1988|) . The specific in- 
tensity, Iuj{t), is calculated using = n p h(s, t)ech/4ir, 
where e — fku. P(u,^/) is the total power emitted by an 
electron having Lorentz factor 7 per unit frequency cj, 
and is given in H2.1.1I 

The time derivative of photon distribution due to syn- 
chrotron emission is given by 



R s (e,t) 



djP(u, 7)n e ±(7, t), 



(G) 



where n e ±(-f,t) = n e -(7,t) +n e +(7,t). 

In a homogeneous plasma, the self absorption coeffi- 
cient is given by 

^2 



a(e, t) 



8m.w 2 



d 1 P(to^)pr 



d_ 



n e ±(7,t) 



Pi 2 



(7) 



llGinzburg fc Svrovatskii I 119691 iRvbicki fc Lightman I 
11979ft . 

2.1.1. Cyclo- Synchrotron emission spectrum 

The power (energy /time/sr/frequency) emitted by a 
single electron moving with velocity (3 = v/c in a fre- 
quency range lo..uj + dto, at an angle 9 with respect to 
the magnetic field, is given by 



V°° ( 00s 0— /3|| \ 2 j2 ( s 
Z^m=l I sinff J • J m\- h ) 



(8) 



S(y)duj 



/see iRekelH 17963 IFTuTzhiirg fc Svrovatskii I MM 
iMahadevan. Naravan fc Yi Ill996j) . Here 



x = — f3± sm9, 

UJ 



(9) 



UJb 

^0 = — ; 

7 



qB_ 

m e c 



y = mujo — lu(1 — (3u cost 



(10) 
(11) 

J m (x) is the Bessel function of order m, J' m {x) its deriva- 
tive, = (3 cos Op and /3j_ = /3sin# p are the velocity 
components parallel and perpendicular to the magnetic 
field, and 6 P is the angle between the electron velocity 
direction and the magnetic field. The presence of a 5 - 
function implies that the emission occurs at discrete fre- 
quencies. 

The total power emitted by a single electron having 
Lorentz factor 7 = (1 — /3 2 )" 1 / 2 per unit frequency uj is 
given by integrating r? u (/3, 9) over the solid angle d£l = 
sin 9d9d<j). For an isotropic distribution of electrons, the 
mean radiated power is given by 

^,7) = Jf; = J^#p/o<co S p ) 

x / o 7r d0/_ 1 (i(cos6')7? w (/3, 9) 
= 2?r d(cos9 p ) J^ 1 d(cos9)i] u ((3,9), 

(12) 
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where the factor of l/4ir comes from angular normaliza- 
tion of the isotropic distribution, and the factor of 2 is 
due to integration on half of the range of cos 9 p . 

In the synchrotron limit 7 > 1, the Bessel functions 
can be approximated by modified Bessel functions, re- 
sulting in the well known result 



V3q 3 B sin 8 P 
27rm P c 2 



F(X), 



where 



X = — , uj c — -7 sin 



3 2 <i B ■ 

-7 SllLVp. 

2 m e c 



F(X) is given by 

F(X) = X [ K 5/3 (OdC, 
J x 



(13) 
(14) 

(15) 



where -^5/3 (£) is modified Bessel function. The function 
F(X) was tabulated in, e.g., iGinzburg fc Svrovatskii I 
(196|). 

The power emitted by a single electron is given by in- 
tegrating eqs. El E| over all frequencies, 



P s ( 7 )ee / P(u J , 1 )doj = 



2g 4 B 2 1 2 (3 2 sin 2 Q p 
3m 2 c 3 



(16) 



(see, e.g.. lRvbicki fc Lightman Ill979|) . 

The calculation method of the cyclo-synchrotron emis- 
sion spectrum P(w,7), is determined by the electron 
energy: (i) For low energy electrons, having 7 < 3.2 
(/3 < 0.95), integration of equation 1721 is carried out ex- 
plicitly, at all frequencies up to u)/u>b = 200. Above this 
frequency, no emission is assumed, (ii) For electrons with 
3.2 < 7 < 10, equation 1121 is solved up to ui/ujb = 100. 
Above this frequency, the approximate synchrotron spec- 
trum (cq. IT3|) is calculated up to w < 10w c . (hi) At high 
electron energies, 7 > 10 the synchrotron spectrum (eq. 
I13|) is calculated in the range 0.001w c < u> < 10ui c . 

2.2. Compton scattering 

The total power emitted by Compton scattering by a 
single electron having Lorentz factor 7 into a unit vol- 
ume, is given by 



Pc(l,t)= J dai J da 



d 2 N( 1 ,ai) 
dtda 



n p h(ai,t)(a - ai), 
(17) 

where d 2 N "(7 , ai) / 'dtda is the rate of scattering by a 
single electron having Lorentz factor 7 passing through 
space filled with a unit density (1 photon per unit vol- 
ume), isotropically distributed, mono-energetic photons 
with energy a\rn e c 2 . Note that the Compton power can 
be negative (i.e., the electron gains energy), depend- 
ing on the initial photon number density distribution, 
n p h(ai,t). 

The time evolution of the photon number density due 
to Compton scattering is given by 



Rc(a,t) 



J d-yn e ±(-f,t) 
x J dain ph {ai, 



t) 



d 2 N(-y.a 1 ) 
dtda 



d 2 N(~f,a) 



dtda 1 



(18) 



2.2.1. Compton scattering spectrum 

The rate of scattering by a single electron having 
Lorentz factor 7 passing through space filled with a unit 
density, isotropically distributed, mono-energeti c pho- 
tons with energy aim e c 2 was first derived by I.Tones I 
(1961, 



d 2 N( 1 ,a 1 ) 



nr^ca 



[f(c + )-f(c_; 



(19) 



dtda 27 4 /3a 2 

Here, a is the energy of the outgoing photon in units of 
m e c 2 , ro is the classical electron radius, f3 = (1 — 1/7 2 ) 1 ' 2 , 
£± are the upper and lower integration limits (see below) 
and F(() is given by the s um of 1 2 fun ctions obtained 



by solving equation (21) of 
Solving equation (21) of 

by 



Jones 



Jones 



1968) 



:S4 



F(() is given 



f(C) = E^i/i. 

fx 

h 
h 
h 

h 
h 



s/Ei 
a( 



ai _ 

7 a 



/EC 




< T \ (a 1 / 7 + l)- 



!io g (^ 



h=(l) 2 1 VE- 2 /c-(l) 2 d 



h 

/10 = 

hi = 

/l2 = 

Here, 



E'2 

c ■ 



ay 



la 



( £LL 
V T 



017 



2 2c 



aiY 



'(V^/d) c>0; 
V-cC/d) c<0, 

sinh -1 c^/d 
(-l)sin- 1 ' 



2C 



|C|3 



</d 



c> 0; 
c < 0, 



sinh 



M 3 



(-l)sin '(^c, 



c> 0; 
c < 0. 



(20) 



a = l/7 2 (ai+7) 



1 



1, 



(21) 



6 = 2ai/7, 
c = (7 — a) 2 - 
<i = 2a/7, 
Pi = a - 6C, 
P 2 = cC 2 + dC- 

The integration limits depend on the energy of the out- 
going photon, a, best presented as a function of the pa- 
rameter p ee a/ ai. The minimum value of p is 5 

1-/3 



Pn 



l + f3 + 2a 1 /~/' 



(22) 



3 Note that in eq. 21 of Jones, there is a misprint by a factor of 
a in the one before last term. 

4 Note, though, that Coppi & Blandford 1 1990) claim about an 
error in eq. (20 ) ofUoiSH lffi}6Sl) is incorrect. In fa ct, eq. Al.l o f 
ICoppi fc Blandford I JT990T) is identical to eq. (20) of IJones I J1968I) . 

5 Note tha t there is a misprint in the result that appears in 
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while the upper value of p is limited by the kinematics, 

ftoax,i = l + (7-l)/ai, (23) 
and by the requirement that £<!+/?, 



Pmax,2 



(24) 

resulting in p max = min(p 

max,l j /^max,2 ) • 

For a given p, p min < p < p max , the integration bound- 
aries are 

C-(p) = maxfp|(l + Q!i/7-/5Q!i/7) 

- [(1 + a x / 7 - pai/7) 2 - V7 2 ] 1/2 },!-/?), 

(25) 

and 

C+0) = mm^|(l + Q!i/7-pQ!i/7) 

+ [(1 + a x / 7 - pai / 7 ) 2 - 1/7 2 ] V2 }, 1 + /?) • 

(26) 

For an energetic electron, 7 ^> 1 and 7 > aj, equation 
[TU1 can be simplified and the scattering r ate is given by 
(Jones 1968: IBlumenthal fc Gouldlli97fj 



d 2 N(~f, ai ) 
dtda 



2glogg+(l + 2g)(l-g) 



1 JAayyqf /, _ \ 
+ 2 (l+4ai 7? ) V 1 1) 



(27) 



where g = a/Aa\"f 2 {l — a/7) is limited to I/47 2 < g < 1. 

Calculation of the spectrum resulting from Compton 
scattering is determined by the electron Lorentz factor 
7 and the incoming photon energy, a\. (i) For 7 > 10 4 
and ot\ < 10~ 5 , the approximate spectrum (eq. I27|) is 
used, (ii) For all other values of 7, oi, the exact spectrum 
(eq. IT9|) is calculated. The results are tabulated in a 3-d 
matrix (initial electron energy x initial photon energy x 
final photon energy) , and are used in calculating the time 
derivatives of electron and photon number densities. 

2.3. Pair production 

The production rate of particles having Lorentz fac- 
tor in the range 7.7 + by an isotropic photon 
field with photon dens i ty n P h(a,t), was calculated by 
IBotcher fc Schlickeiser I Ipfll. 



9w eP (7,t) _ 3 

— at — ~ 4 



C J0 dQ! l a* 



poo 

x / ; 1 



1 , , 1 da 2 pfa V 1 
|^_ s7 +i- ai j ^ otl 



(28) 



ai t 2 are the scattering photons energies in units of m e c 2 , 
E = ct\ + ct2, and a cm is the photons energy in the 
center of momentum frame, given by 2a 2 . m = ct\ ■ a 2 . 



The functions H± are calculated using 
c± = ("1,2 - 7) 2 - l , 
z a\ 2 + aia 2 ± 7(02 - a\) 



(29) 
(30) 



For c± 0, H± are given by 



8^ 



CtiCt2 



c ± a cm 



1 + Ql/7 + y/(l + ai/7) 2 ~ 1 + (3 2 - 2ai/7(l + P) 

1-/3 + 207/7 



Q; lQ ; 2 — 1 



ai«2 
4 



2 



Ct\OL-2. 



&crn 



(31) 




In facm-y/Ci + V a l a 2 + C±( 

7% arcsin («<™^-^) 



(aia 2 ) 3 / 2 
1_ 

ta c „ 



v /Q i Q 2 ' 



The upper and lower integration limits a 
given by 



u 

era i 



c± > 0, 

c± < 0. 

(32) 

(33) 
are 



where 



max 



(34) 



, a, 6 \ 



i«c™J = 2 - 7] + 1 ± V(tP - 7] + 1) 2 - # 2 



(35) 

The total loss rate of photons in the energy range 
a\..a\ + da\ by pair production is given by 



R P (a\,t) = -n ph (oti,t)j Jd(cos0)(l 



COS( 



1: 



da 2 n p h(a 2 , t)a(a\, a%, 0), 



where 



and 



(36) 



tWW 2 ) 

2/3'(/3' 2 -2) + (3-/3' 4 )ln(i±|; 







1 



1/2 



(37) 



(38) 



0:10:2(1 — cos 6 

(ICould fc Schreder I Il967t lLanTI Il999t) . The result- 
ing particle spectra are symmetric for electrons and 
positrons. 

Calculation of the photon loss rate is carried out us- 
ing equation 1361 The spectra of the emergent pairs is 
calculated in accordance to the photons energies: (i) For 
1.001 < 0:1-0:2 < 10 4 , eauationl28lis solved, and the exact 
spectrum is obtained, (ii) For a.\ ■ a 2 < 1.001, monoen- 
ergetic spectrum of the created particles assumed, with 
energy (ai + a 2 )/2. (iii) If 0:1 • o 2 > 10 4 , one of the cre- 
ated particles energy is taken to be o max = max(ai, 02), 
and for the second particle the energy is approximated 
as o min + l/(2o min ), where a min = min(ai,a<2). 

2.4. Pair Annihilation 

The total loss rate of electrons having Lorentz factor 
71.. 71 +e?7i due to pair production (in the plasma frame) 
is given by 



dn 



,(7i,*) 



dt 



-( 7 i,t)3^/d(coBfl) 



I dn , 



(39) 
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where 7 2 = 7172(1 + P1P2 cos 9) is the positron Lorentz 
factor in the electrons rest frame, (3' 2 is its velocity in this 
frame and dn' /dn = 71(1 + /3i/?2 cos 9). The cross section 
for a positron having Lorentz factor 7 to annihilate with 
an electron at rest, <J a nn{l), is given by 



t (7) 



3 or 
8 7+1 



X 



7+2 



(40) 



ijSvensson 1119821 ILang l l 1999) . The loss rate of positrons 
is calculated in a similar way. 

The annihilation rate is ca lculated by solvin g equation 
1391 Since (i) It was shown in ISvensson I ((19825) that for a 
large region of 71, 72 the photons spectrum is narrowly 
peaked around £1.2 = 71.2'TieC 2 , and (ii) We found nu- 
merically that calculation of the exact particle spectrum 
resulting after pair production, compared to the approx- 
imate particles spectrum 71,2 = Ei^/rrieC 2 , did not have 
a significant effect on the resulting photon spectra, we 
decided not to include calculation of the pair annihilated 
photon spectra in this version of the code. The emergent 
photons energies assumed to be equal to the reacting 
particles energies, £1,2 = 7i,2«ieC 2 , thus 



-ym e c 2 ,t) 



dn eA (j,t) 
dt 



(41) 



2.5. Photon and positron production by it decay 

Photo-meson interactions between energetic protons 
and low energy photons result in production of 7r's. The 
fractional energy loss rate of a proton with Lorentz factor 
7 P due to pion production is 

L TT KIP' 1 ) — 7p dt 

2r^ c JT dea 7T (e)((e)e J™ 2 ^ dxx~ 2 n ph (x,t) 



(42) 

ijWaxman Ez Bahcall Ill997j) where ov (e) is the cross sec- 
tion for pion production for a photon with energy e in 
the proton rest frame, £(e) is the average fraction of en- 
ergy lost to the pion, and £0 = 0.15 GeV is the threshold 
energy. For a flat photon spectrum {e 2 n p h{s) ex e a with 
a ~ 0), the contribution to the first integral of equation 
l4*2l from photons at the A-resonance is comparable to 
that of photons of higher energy, thus 



dxx 2 n p h{x,t), 



: 5 x 10 28 cm 2 and £ peak c 
= 0.3 GeV , and Ae 



(43) 
0.2 at the 
; 0.2 is the 



where a peak 
resonance e 
peak width. 

The rate of proton energy transfer to pions is given by 

P„{j,t) = t-\ lp ,t) x lp m p c 2 . (44) 

The energy loss rate of protons is calculated by numer- 
ical integration of the integral in equation 031 This cal- 
culation is carried out only in those cases where the A- 
resonance approximation is valid, and can easily be ex- 
tended to any photon spectrum by explicit integration of 
the integrals in equation 021 Roughly half of this energy 
is converted into high energy photons through the tt° de- 
cay. Each of the created photons carry 10% of the initial 



proton energy, thus the photon production rate is given 

by 

Rtt(s = l P m p c 2 1 '10, t) = 5t^ 1 (7 P ,t)n p (7p,i), (45) 

where n p (^ p ) is the number density of protons at energy 
jpm p c 2 . Half of the energy lost by protons is converted 
into 7r + , that decays into positron and neutrinos, n + — > 
/i + + 17, — > e + + v e + + v^. The 7r + 's energy is roughly 
evenly distributed between the decay products, thus the 
positron carries 5% of the initial proton energy, and the 
positron production rate is given by 

QAl = 7 P K>/™ e )/20, t) = 2M-\ lpi t)n p ( lp , t). (46) 

Equations 1451 1461 provide only a crude approximation 
to the spectrum of high energy photons and positrons 
produced by pion decay. However, photons and positrons 
that are created by pion decay are typically very ener- 
getic, and participate in the high energy electro-magnetic 
cascade. Since these particles and photons' energy is 
spread among the cascade products, and the final cas- 
cade spectrum has only weak dependence on the initial 
spectrum, it is appropriate to use the approximate ex- 
pressions in equations 1451 and 1461 

3. NUMERICAL APPROACH 

Several integration methods are used in solving the ki- 
netic equations. Simple, first order difference scheme was 
found adequate, except when dealing with synchrotron 
self absorption and with the evolution of the rapid high 
energy electro-magnetic cascades. Synchrotron self ab- 
sorption calculations are carried usin g Cranck-Nickolsq n 
second order integration scheme (see lPress et ]Ol992). 

The particle and photon distributions are discretized, 
spanning the energy range relevant to the problem. Note 
that in the problems involved, this energy range can ex- 
tend over 20 decades (see 21 below) . Spectra of cyclo- 
synchrotron emission (eqs. El Compton scatter- 
ing (eqs. El E3 and pair production (eq. EHJ are pre- 
calculated and stored in tables. 

Following simultaneously the evolution of the rapid 
high energy electro-magnetic cascade and the much 
slower evolution of low energy processes, is difficult. The 
'stationary' approximation used in previous works in 
treating the evolution of high energy particles (see, e.g., 
Fabian et al. 1119861 lllghtman fc Zdziarski 119871 iCoppi I 
1992) can not be used, due to the non-linear nature of the 
cascade: As an energetic particle loses its energy, many 
secondaries are created, which, in turn, serve as primaries 
for further development of the cascade. As the cascade 
evolution is powered by inverse Compton scattering, pair 
production and annihilation, the injection rate of ener- 
getic photons and pairs depends on the entire particle 
and photon spectra. 

Therefore, in treating this problem, a fixed time step is 
chosen, typically 10 -4 5 times the dynamical time. Nu- 
merical integration is carried out with this fixed time 
step. At each time step, the cascade evolution is followed 
directly. Direct numerical integration is carried out only 
for the electrons, positrons and photons for which the 
energy loss time or annihilation time are longer than 
the fixed time step. Electrons, positrons and photons 
for which the energy loss time or annihilation time are 
shorter than the fixed time step, are assumed to lose all 
their energy in a single time step, producing secondaries. 
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The secondaries' spectra are determined by the spectra 
of the various physical processes, as presented in and 
by the relative rates of these processes. We discriminate 
between high energy secondaries, which are secondaries 
for which the energy loss time or annihilation time are 
shorter than the fixed time step, and low energy secon- 
daries which lose their energy or annihilate on a time 
scale longer than the fixed time step. The calculation 
is repeated for the high energy secondaries, which are 
treated as a source of lower energy particles, until all 
the cascade energy is transferred to low energy particles. 
Since in each step of the cascade calculation, part of the 
energy is transferred into low energy particles which do 
not participate in the cascade, convergence is guaran- 
teed. In order to check for convergence of this method, 
we repeat the complete calculation with a shorter time 
step. 

The time derivative of particle distributions due to 
synchrotron emission and Compton scattering is calcu- 
lated by solving the continuity equation, dn(pi,t)/dt = 
dj(-f,t)/dE, where j(j,t) = n(j,t)P(j), and P(7) is the 
emitted power (see eq. 0. In solving this equation, flux 
limiter is used to ensure convergence for large time steps, 
and Neumann boundary conditions for the flux t) at 
the boundary points are used. The rate of change of 
particle distributions due to pair production and annihi- 
lation are calculated using eqs. EHl QUI Conservation of 
particle number and energy is forced using Lagrange mul- 
tiplier method. This method was found to allow faster 
convergence (larger time steps). 

At the low end of the particle spectrum, electrons and 
positrons gain energy via direct Compton scattering, on 
a timescale shorter than the fixed time step. In parallel 
to gaining energy, these particles lose energy via syn- 
chrotron emission on a much longer time scale, thus pro- 
viding another challenge to numerical integration. Defin- 
ing 'very low energy particles' as particles that gain en- 
ergy on a time scale shorter than the fixed time step, we 
treat this problem in the following way. At each time 
step, calculation of the number density of these parti- 
cle is repeated iteratively, until the particle distribution 
converges, and the total emissivity equals the absorption. 
At each of the iteration steps, the calculated emissivity 
and absorption are stored, and used in the calculation of 
the photon emission from these particles. Convergence 
of this method as well is checked by repeating the calcu- 
lation with smaller time step. 

4. EXAMPLES OF NUMERICAL RESULTS 

We give below several examples of the results of numer- 
ical calculations of GRB prompt emission spectra. De- 
tailed description of numerical results of prompt emission 
spectra and early afterglow emission spectra are found in 
iPe'er fc Waxmanl (2004b c). Our calculations ar e done 
in the framework of the fireball mode l (see, e.g.. iPiranl 
120001 iMeszaros 1^0021 iWaxman ll2003|) . where the emis- 
sion results from electron acceleration to ultra-relativistic 
energies by internal shocks within an expanding wind. 

4.1. Basic assumptions, plasma conditions and particle 
acceleration 

We calculate the emergent spectra following a single 
collision between two plasma shells. Denoting by F the 
characteristic wind Lorentz factor, and assuming varia- 



tion Ar/r ~ 1 on a time scale At, two shells collide at 
radius ri = 2T 2 cAt. For Ar/r ~ 1, two mildly relativis- 
tic (r s — 1 ~ 1 in the wind frame) shocks are formed, 
one propagating forward into the slower shell ahead, and 
one propagating backward (in the wind frame) into the 
faster shell behind. The comoving shell width, measured 
in the shell rest frame, is Ai? = FcAt, and the comoving 
dynamical time, the characteristic time for shock cross- 
ing and shell expansion measured in the shell rest frame, 
is tdyn — TAi. The shock waves, which propagate at 
relativistic velocity v s ~ c in the plasma rest frame, dis- 
sipate the plasma kinetic energy and accelerate particles 
to high energies. Since the shock velocity is time inde- 
pendent during td yn , the shock-heated comoving plasma 
volume is assumed to increase linearly with time, i.e., 
constant particle number density is assumed. 

Under these assumptions, the shocked plasma condi- 
tions are determined by six free parameters. Three are 
related to the underlying source: the total luminosity 
L = 10 52 L52 ergs -1 , the Lorentz factor of the shocked 
plasma, T — 10 2 5 r2.5, and the variability time At 
1 Ai_4 s. Three additional parameters are related 
to the collisionless-shock microphysics: The fraction of 
post shock thermal energy carried by electrons e e = 
10~ a5 e e! o.5 and by magnetic field, es = 10~ - 5 e_B^o.5, and 
the power law index of the accelerated electrons' Lorentz 
factor distribution, d log n e jd log 7 = — p assumed to ex- 
tend over the range 7 min < 7 < 7 max . 

The comoving proton number density is 

L 



4irr 2 cT 2 m v c 2 



= 6.7 x io 14 l 52 At: 



cm 



-3 



The internal energy density is Ui n t 
suiting in a magnetic field, 



(47) 

L/(47rr 2 cL 2 ), re- 



B = 



2L 6 c 2 At 2 



= 2.9 x 10 6 LH 2 e 1/2 



-1-3 

B,-0.5 x 2.5 



r 2 -j Atzi G. 

(48) 



4.1.1. Particle acceleration 



Since the details of the acceleration mechanism are 
not yet known, we adopt the common assumption of 
a power law energy distribution of the accelerated elec- 
trons Lorentz factor, 7. The power law index p of the 
accelerated particles is a free parameter of the model. 
The maximum Lorentz factor of the accelerated elec- 
trons, 7max is obtained by equating the acceleration time, 
tacc — ^m e c 2 1 1 cqB , and the synchrotron cooling time, 



syn 



9m^c 5 /4g 4 B 2 7, to obtain j n 



{Qitq/ o~tB) 



1/2 _ 

The accelerated particles assume a power law energy dis- 
tribution above a minimum Lorentz factor 7 m i n , which is 
obtained by simultaneously solving 



(49) 



where n e and u e = e e Ui n t are the number and energy 
densities of the electrons. The injected particle distri- 
bution below 7 m j n is assumed thermal with temperature 
9 = kT/m e c 2 — 37 m in, and exponential cutoff is assumed 
above 7 max . In the results shown below, no proton ac- 
celeration is assumed. 

Our calculations are carried in the plasma (comov- 
ing) frame. The particle distributions are discretized, 
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spanning total of 10 decades of energy, (7/3 m i n = 10 3 
to 7/? max = 10 7 ). The photon bins span 14 decades 
of energy, from a min = e min /m e c 2 = 10~ 8 to a max = 
£max/™ e c 2 = 10 6 . No a-priori photon field is assumed. 

4.2. Low compactness 

We examined the dependence of the emergent spec- 
trum on the uncertain values of the free parameters of 
the model. We found that the spectral shape strongly 
depends on the dimensionless compactness parameter I, 
defined by I = LaT/Rm e c 3 , where L is the luminos- 
ity, and R is a characteristic length of the object. For 
low value of the comoving compactness, I' < 10, the 
optical depth to pair p roduction and to scatter ing by 
pairs is smaller than 1 ijPe'er fc Waxman"! 120 04b). thus 
synchrotron-SSC emission model provides fairly good ap- 
proximation of the resulting spectrum. Therefore, before 
applying our model to examine realistic scenario, (i.e., 
comparison with observations), we first compare our nu- 
merical results to the analytical model predictions, in 
parameter space region where the later are valid. Figure 
^presents numerical results in this case, where a power 
low index p = 3 was used to allow the synchrotron and 
Compton peaks to be distinctively apparent. 

The synchrotron peak presented in Figure^at £° b eak — 
10 keV is in excellent agreement with the analytical re- 
sults of the optically thin synchrotron model prediction, 

ob. = h 3_gB_ 2 r 

c peafe '"2 m e c 'mm 1 

= 1.4 x io 4 l% 2 < _ o 5 e V2_ o 5 r -g At -i eV; 

(50) 

where 7 m i n ~ e e (m p /m e )(p — 2)/(p — 1) was used. The 
Lorentz factor of electrons that cool on a time scale that 
is equal to the dynamical time scale is j c ~ 1, thus above 
£ V eak the spectral index, vF v cx v a , is a = 1 — p/2 = 
— 1/2, while below e pea k, a = 1/2. The self absorption 
frequency, e°^ a ~ 100 eV , is somewhat lower than the 
self absorption frequency predicted for a pure power law 
distribution of the electrons, 

ef sa = 600 L% 3 e\%_ 5 e^_ . 5 T 2 J /3 AC* eV, (51) 

where a power law index p — 2 for particles below 7 m i n 
is assumed. This discrepancy is due to the fact that 
the low energy particles are not power-law distributed, 
but have a quasi-Maxwellian distribution due to photon 
reabsorption (see figure 0). 

Without pair production, the SSC model predic- 
tions of the Compton scattering peak, at SjQ peak = 

7min e peafe = 1-3 GeV , agrees well with the numerical re- 
sult efc peak = 1.5 GeV . The 1 GeV flux is comparable 
to the flux at 10 keV, as predicted by analytic calcula- 
tions based on the Compton y parameter, y = 1 in the 
scenario presented in figure ^ 

Pair production causes a cutoff at high energies. For 
a flat spectrum, e 2 n p h{£) oc e° (which is a good approxi- 
mation provided es is not much below equipartition), the 
optical depth to pair production is well approximated by 

r 77 (e) = ARn ph (s)-^a T 

lo g ( ) (nieC 2 ) 2 16 1 ' 

V £ peak J 



and is larger than unity at 

if > 3 x 10 8 log ( e f^) L^c-io.^' 6 Ai_ a eV, 

(53) 

in an excellent agreement with the numerical results. 
Here U p h is the photon energy density, given by U p h ~ 
e e i/47rr 2 cr 2 . For this value of the compactness, pair 
annihilation does not play a significant role, while scat- 
tering by the created pairs flattens the spectrum at 
10 keV - 1 GeV . 

Even though the analytic approximation is in a fairly 
good agreement with the numerical calculations, there 
are important discrepancies between the analytic approx- 
imation and the numerical calculation. The electron dis- 
tribution shows a peak at 7 ~ 1.05 ((3 ~ 0.3), result- 
ing in a deviation of the self absorption frequency from 
the analytic calculation. These electrons affect the high 
energy spectrum by Compton scattering, resulting in a 
nearly flat {a ~ 0) spectru m above 10 keV . We showed 
l)Pe'er fc Waxmanl 12004111 that the spectrum is nearly 
independent on the power law index of the accelerated 
electrons, p. 

4.3. High compactness 

Figure |3 shows an example of our numerical results for 
large comoving compactness, V = 250. At large value 
of the compactness parameter, V > 30, the synchrotron- 
SSC model predictions do not provide an appropriate 
description of the spectrum. Therefore, the numerical 
results may provide some insight on the inconsistency 
between some of the observations and the analytical pre- 
dictions, as mentioned in ^ 

In the scenario of large compactness, Compton scat- 
tering by pairs becomes the dominant emission mech- 
anism. Both electrons and positrons lose their energy 
much faster than the dynamical time, and a quasi- 
Maxwellian distribution with an effective temperature 
8 = kT/m e c 2 ~ 0.05 — 0.1 is formed. Photons upscat- 
tered by the pairs create the peak at T8m e c 2 ~ 5 MeV . 
The results shown in Figure are not corrected for the 
fact that the optical depth to scattering by pairs is large, 
t± ~ 10 (see jJSJ. Therefore, the emergent spectral 
peak is expect to be at lower energy, at ~ 1 Me V (for 
detailed discussion sec IPe'er fc Waxlnanl I2004bl) . The 
moderate Compton y parameter, y ~ A9t w \B^\T\ re- 
sults in a spectral slope vF v oc v a with a s» 0.5 be- 
tween e ssa w 3 keV and e pea k ~ 5 MeV . The peak at 
Tm e c 2 ~ 10 2 r2.5 MeV is formed by pair annihilation. 
The self absorption frequency, e ssa — 3 keV is well be- 
low the prediction for a power law index p = 2 of particles 
below 7mi n . This is attributed to the quasi-Maxwellian 
distribution of particles at low energies (see Figure |2J ■ 
The inverse Compton peak flux is lower than the syn- 
chrotron peak flux, due to the Klein-Nishina suppression 
at high energies. 

Even though the results presented here are for illus- 
trative purposes only, and are not aimed to explain a 
particular observation, we note that the obtained nu- 
merical results are is agreement with some of the observa- 
tions, that were found inconsistent with the optically thin 
synchrotron- SSC model predictions . Examples are the 
steep slopes obs e rved a t low energies llPreece et al. Il998t 
iFrontera et~aT\ 120001: iGhirlanda. Celotti fc Ghisellini I 
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FlG. 1. — Time averaged GRB prompt emission spectra obtained 
after 2 shells collision, characterized by low compactness param- 
eter. Results are shown for L = 10 5 erg, e e = eg = 10~ 5 , 
p = 3, At = 10~ 2 s , r = 300. The comoving compactness pa- 
rameter is I' = 2.5. Dotted curve: cyclo- synchrotron emission 
only. Dash dotted curve: synchrotron emission and self absorption 
only. Dashed: synchrotron emission, synchrotron self absorption 
and Compton scattering. Solid: All processes included, including 
pair production and annihilation, but excluding proton accelera- 
tion. Luminosity distance d^ = 2 X 10 28 and z = 1 were assumed. 




10° 10 1 10' 10" 
electrons momenta: yx p 

FlG. 2. — Particle distribution at the end of the dynamical time. 
Thick: At = 10" 4 s, I' = 250; Thin: At = 10" 2 s, /' = 2.5. 
All other parameters are the same as in Figure Solid: electron 
distribution, dash-dotted: positron distribution. The dotted line 



shows a Maxwellian distribution at temperature 9 
0.08. 



kT/meC 2 



12003^ . and the steep slopes above s° b eak obtained by 



iBaring fc Brabv I (I2004T) . Further results o f our study 
are presented in iPe'er fc Waxmanl l)2004bf) . Compari- 
son of the numeric al results with the hig h energy com- 
p onent reported bvlGonzalez et al. I i|2003|) are presented 
in IPe'er fc Waxman I l)2004aj) . 



5. SUMMARY & DISCUSSIONS 

We described a time dependent numerical model which 
calculates emission of radiation from relativistic plasma, 
composed of homogeneous and isotropic distributions of 
electrons, positrons, protons and photons, and perme- 
ated by a time independent magnetic field. We assume 
the existence of a dissipation process, which produces 
energetic particles at constant rates. The particles in- 
teract via cyclo-synchrotron emission, synchrotron self 
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FlG. 3. — Time averaged GRB prompt emission spectra obtained 
after 2 shells collision, characterized by high compactness param- 
eter (not corrected for the high optical depth to Thompson scat- 
tering). Results are shown for L = 10 52 erg, e e = eg = 10 -0 5 , 
p = 3, At = 10 -4 s , r = 300. The comoving compactness pa- 
rameter is I' = 250. Dotted curve: cyclo- synchrotron emission 
only. Dash dotted curve: synchrotron emission and self absorption 
only. Dashed: synchrotron emission, synchrotron self absorption 
and Compton scattering. Solid: All processes included, including 
pair production and annihilation, but excluding proton accelera- 
tion. Luminosity distance = 2 X 10 28 and 2 = 1 were assumed. 



absorption, inverse and direct Compton scattering, e ± 
pair production and annihilation, and photo-meson inter- 
actions which produce energetic photons and positrons 
following the decay of energetic pions. Exact cross sec- 
tions valid at all energies, including the Klein-Nishina 
suppression at high energies, are used in describing the 
physical processes. Exact spectra are used in describing 
cyclo-synchrotron emission, synchrotron self absorption, 
Compton scattering and pair production, and approxi- 
mate spectra are used in the description of pair annihi- 
lation. 

We explained in our unique integration method, 
which overcomes the challenge of the many orders of 
magnitude difference in characteristic time scales. We 
presented the various integration techniques used for 
solving the kinetic equations describing the evolution of 
particle and photon distributions at all energy scales. By 
following directly the development of the rapid, high en- 
ergy electro-magnetic cascades at each time step, we ob- 
tain the spectrum at high energies, up to > 100 TeV . 
Our method enables to follow the development of the 
spectrum created in the parameter space region of large 
compactness, where no analytic approximation is valid. 
This method also improves over previous ones by provid- 
ing a more accurate treatment of photon emission and 
absorption in the presence of magnetic fields. 

We have given several examples of numerical calcula- 
tions in 21 In parameter space regions where analytical 
approximations are valid, our numerical results are in 
good agreement with analytic results. We have pointed 
out some significant discrepancies between the analyt- 
ical approximations and the numerical calculation, and 
explained their origin. We presented examples of new re- 
sults for parameter space regions where analytic approx- 
imations are not valid. We pointed out that our results 
are consistent with numerous observations, including ob- 
servations that are inconsistent with the optically thin 
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synchrotron-SSC model predictions. Further results of 
our study of GRB prom pt emission and early after glow 
emission can be found in lPe'er fc WaxmaiTI l)2004bl lcl) . 

The next generation high energy detectors, such 
as SWIFT 6 and GLAST 7 satellites, and the sub-TeV 
ground based Cerenkov detectors, such as MAGIC 8 , 
HESS 9 VERITAS 10 and CANGAROO III 11 are expected 
to increase the GRB prompt emission and early afterglow 

6 http://www.swift.psu.edu 

7 http://www-glast.stanford.edu 

8 http:/ /hegral. mppmu.mpg.de/MAGICWeb 

9 http://www.mpi-hd.mpg.de/hfm/HESS/HESS.html 

10 http:/ /veritas. sao.arizona.edu/ 

11 http:/ /icrhp9. icrr.u-tokyo.ac.jp/ 



emission detection rate by an order of magnitude, to al- 
low detection of > GeV emission from GRB's, and to 
detect the high energy spectra of thousands of AGN's at 
various distances. Thus, detailed numerical models, that 
are capable of producing accurate spectra over a wide en- 
ergy scale, are necessary for analyzing and understanding 
the experimental data. 



APPENDIX 
COULOMB SCATTERING 

In the limit of relativistic particle scattering off cool thermal pair distribution (0< 1,7> 1), the energy loss rate of 
the relativisti c particle can be approximated by d'y/dt ~ — 3/2axcn± In A, or t~J- = —{l/"f)d^/dt ~ 4ncr 2 l n±j~ 1 In A 
ijGould 11197,^ - where n± is the number density of the thermal pairs. The relevant value of the Coulomb logarithm A, 
is A « 7 1 / 2 m e c 2 //io;p, where uj p — (47rn±e 2 /m e ) 1 / 2 is the plasma frequency. 

Assuming that the pairs' energy distribution is thermal, that their number density is n± = fn p where n p is the 
proton number density, and that n p rj Ui n t/m p c 2 where Ui n t is the internal energy density, comparing the Coulomb 
cooling time and the synchrotron cooling time, t syn = Qm^c 5 / Aq A B 2 ~{ , using B 2 = 8TreBUi n t, gives 

t syn _ 9 m e / In A 1 f -l f 41 \ 
t ee ~ 8m p e Bl 2 ~ 3 7 2 h "^-o^ ^) 

where typical values / — 10/i, and In A « 20 assumed. It is therefore concluded that for relativistic electrons, and for 
magnetic field not many orders of magnitude below equipartition, electrons lose their energy by synchrotron emission 
on a time scale much shorter than the e nergy loss time by Coulomb scattering. A more accurate approximation of t ee 
(Haug 1988; Coppi & Blandford1ll990D . does not change this result. 
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